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Abstract 

The stability of spin-glass (SG) phase is analyzed in detail for a fermionic 
Ising SG (FISG) model in the presence of a magnetic transverse field V. 
The fermionic path integral formalism, replica method and static approach 
have been used to obtain the thermodynamic potential within one step 
replica symmetry breaking ansatz. The replica symmetry (RS) results 
show that the SG phase is always unstable against the replicon. More- 
over, the two other eigenvalues A± of the Elessian matrix (related to the 
diagonal elements of the replica matrix) can indicate an additional insta- 
bility to the SG phase, which enhances when V is increased. Therefore, 
this result suggests that the study of the replicon can not be enough to 
guarantee the RS stability in the present quantum FISG model, espe- 
cially near the quantum critical point. In particular, the FISG model 
allows changing the occupation number of sites, so one can get a first 
order transition when the chemical potential exceeds a certain value. In 
this region, the replicon and the A± indicate instability problems for the 
SG solution close to all range of first order boundary. 

1 Introduction 

The interplay between disorder, frustration and quantum effects presents several 
challenging issues. For instance, the infinite-range Ising spin-glass model in a 
transverse magnetic field has been studied with various techniques that show 
controversial results. It is well known that the classical Sherrington-Kirkpatrick 
(SK) model [T| presents a continuous phase transition in which the spin-glass 
(SG) phase has the free energy landscape composed by many almost degenerated 
thermodynamic states separated by infinitely high barriers [5] . The controversy 
is whether or not the quantum tunneling between free energy barriers activated 
by the transverse field is able to restore the replica symmetry (RS) in the infinite- 
range Ising spin-glass model. [31 [H [51 O [71 [5] . To answer this question, the 
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usual procedure is to investigate the behavior of the so-cahed repUcon (the 
transversal eigenvalue of the Hessian matrix) to check the local stability of RS 
solution. No much attention is given to the longitudinal eigenvalues since, as 
in the classical case, they would not give any additional important information 
in such an analysis [31 [5]. Nevertheless, there is, at least, one formulation of 
the infinite-range Ising spin-glass problem, the fermionic Ising spin-glass (FISG) 
model [9l [TOl [11] , in which the role of the quantum effects on the longitudinal 
eigenvalues could also be important. 

The FISG model is defined by spin operators which are represented by bilin- 
ear combinations of fermion operators of creation and destruction. That makes 
this model a quite adequate framework to study the competition between spin- 
glass and, for instance, Kondo effect or superconductivity. [TU [131 [14] . Those 
spin operators act on Fock space with four states per site, two of them non- 
magnetic. In fact, this model can also be formulated in two versions [151 [16]. 
In the first one (25' model) , it is used a restriction in the number operator that 
eliminates the contribution of nonmagnetic states. That would be equivalent to 
study the problem in the spin space. In the second version (45 model), mag- 
netic and nonmagnetic states are admitted. In this case, an important aspect 
is the connection between spin and charge correlations [T7]. That can be seen 
clearly, for instance, in the relationship between n (the average of occupation of 
fermions per site) and the diagonal replica matrix elements [131 HZ]. Therefore, 
in the 45 model, variations of charge occupation can influence the magnetic 
properties not only leading the onset of the spin-glass phase, but also changing 
the nature of the phase boundary. For instance, the phase diagram of the 45 
model presents a tricritical point [17j and a reentrance in the first order bound- 
ary phase for specific values of the chemical potential /i [14]. The mentioned 
connection can also appear in the behavior of longitudinal eigenvalues of the 
Hessian matrix. Indeed, in the 45 model these eigenvalues present a non-trivial 
behavior, in which they become complex [131 [17] . 

The presence of a transverse field F produces important changes in the phase 
diagram of 25 and 45 models. For instance, the increase of F leads the freez- 
ing temperature (T/) of both models to a Quantum Critical Point [QCP] at 
F = Fc [Hllll]- Particularly, for the 45-model, when ^ 7^ 0, the increase of F 
not only depresses the second order part of T/ and the position of the tricritical 
point but also destroys the reentrance in the first order part of T/ [T^ . In other 
words, F induces quantum spin flipping mechanisms that affect the replica ma- 
trix elements strongly, in particular, the diagonal ones. Because of that, F could 
interfere in the replicon {Xat) and, it should be remarked, in the longitudinal 
eigenvalues (A±) of Hessian matrix for both 25 and 45 models. 

The previous discussion leads to the question: how exactly would F affect 
A± in these two models? In the 25 model, the diagonal replica matrix elements 
become dependent on temperature when F > 0, which could reflect on the 
stability of the SG phase. Thus, it can be a mechanism to change the behavior 
of A± since, at F = 0, this model reproduces the basic features of the classical 
SK one. One can also remark that F in the 25 model produces effects which 
are limited to the magnetic properties. Nevertheless, that is not the case for the 
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45* one. In this model, there is a mixing between spin and charge correlation 
functions which is important to determine not only magnetic properties [14j 
but also charge ones [13] ■ In this particular model, even for F = 0, the diagonal 
replica matrix elements have a relevant role for the SG solution when fi is 
increased, which can lead to a non-trivial behavior of A± , as discussed previously. 
Consequently, one can expect that the combined effects of variations of T and jjL 
could enhance the non-trivial behavior of \± . Actually, the variations of these 
independent parameters create a complicated interplay between charge and spin 
degrees of freedom in which both magnetic properties and charge distribution, 
for instance, are affected. However, such redistribuition of charge can influence 
the magnetic correlation functions themselves. As a result, the behavior of \± 
in the 45* model can be considerably more complicated than the 25 one. 

Therefore, the main purpose of this paper is to study in detail the RS stability 
for both the 25" and 45" models in the presence of a transverse magnetic field F. 
In particular, for the 45 model, special attention is given to investigate effects 
of combined variations from F and ji on Xat and X±. Furthermore, it is also 
obtained the behavior of n inside the SG phase to clarify how spin and charge 
correlation functions influence each other. In order to accomplish these goals, 
the partition function is obtained in the path integral formalism in which the 
spin operators are given as bilinear combinations of Grassmann fields. The 
thermodynamic potential is found within the static approach {SA) and the one 
step replica symmetry breaking (IS-RSB) scheme, which allows to derive n. 
In particular, to obtain the 25 model, it is introduced the same restraint used 
in Refs. [TOt [TS] to avoid the two nonmagnetic eigenstates. Finally, the local 
stability of the RS solution is studied by obtaining Xat and A± . 

The present work the static approach which neglects time fluctuations of 
spin-spin correlation functions is used [TS]. It is clear that the SA is unable to 
provide reliable quantitative results for the thermodynamics at very low tem- 
perature [19j. However, it has also been shown by Grempel and Rozemberg [20] 
that static ansatz gives reliable results when the temperature is not too low. In 
fact, the results found in Ref. [5D] for the spin-spin correlation function Q{t) 
indicate a behavior which is approximate to the same of their classical counter- 
part within an interval of temperature around a continuous SG transition. On 
the other hand, when a first order transition is present, it would be expected 
that SA could produce naturally reasonable results around the transition |19j . 
This interval of reliability for both kind of transition justifies the use of SA for 
purposes of this work. 

It is important to mention that the scenario described for the 45-model 
when F = has been already found in another model. The classical Ghatak- 
Sherrington (GS) model [U] is a particularly interesting example which displays 
a spin-glass first order phase transition [551|531|211 US] . In the GS model, there 
are three distinct eigenvalues of the Hessian matrix in the limit of number of 
replicas n — > 0. As in the classical SK model the replicon (Xat) is negative 
in the whole SG region. However, other two eigenvalues (A±) can be complex 
or negative [23l [24l [25]. Actually, there is a close relationship between GS 
and FISG models ^7\. In particular, the thermodynamics described in the GS 
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model can also be found in the 4S'-model, such as the reentrance in the first 
order boundary phase. This similarity of the thermodynamics between these two 
models is obtained by a mapping of their spin-glass order parameters equations, 
within the SA and RS solution, given by a relationship between the chemical 
potential fi and the anisotropic constant D of the GS model [13]. In that sense, 
the GS model can be a reference to check the correctness of results obtained in 
the 45"- model when F ^ 0. 

The organization of this article is as follows: in section [51 the FISG model in 
a r field is presented. The thermodynamic potential is obtained within the SA 
ansatz and IS-RSB solution. In this section, it is also done a detailed study of 
the local stability of the RS solution for the present FISG model. Section [31 is 
devoted to describe the numerical results. Section [4| is left for conclusions. 



2 General Formulation 

The infinite-range Ising spin-glass model in the presence of a transverse magnetic 
field is described by: 

H^-J2j.,sts^-2rJ2s^ (1) 

where the sums are run over the N sites of a lattice. The exchange interaction 
Jij among all pairs of spins is assumed to be a random variable with a Gaussian 

distribution 

P(J„) = V^^/(327rJ2) exp(-4.Ar/32 J2). (2) 

The spin operators are defined as [15] : 

S^ = ^[n^^ - n^i] , S'f = ^[4](^il + 4iCt\] (3) 

where fii^ = cl^Ci^ is the number operator, (cia) is fermion creation (de- 
struction) operator and a —'\ or | indicate the spin projections. 

The Hamiltonian given in Eq. ([1]) is defined on the Fock space where there 
are four states per site: one state with no fermion, two states with a single 
fermion and one state with two fermions. Consequently, there are two nonmag- 
netics states. This work considers two models: the 45 model that allows the 
four possible states per site and the 25" model, which restricts the spin operators 
to act on a space where the nonmagnetic states are forbidden. Therefore, the 2S 
model requires a restriction to remove the contribution of these nonmagnetic 
states. It can be obtained by computing only sites occupied by one fermion 
(n^i -|- riii = 1 at every site) in the partition function trace [10[ I15j . 

In this fermionic problem the partition function is expressed by use of the 
Lagrangian path integral formalism in terms of anticommuting Grassmann fields 
{(/) and (/)*) [T5]. The restriction in the 2S-model is obtained by using the Kro- 
necker delta function [Tni[Ii|. Therefore, adopting an integral representation 
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for this delta function, one can express the partition function for both models 
in a compact form: 

Z{y} = e^^''^ J D{(p*(p) \[ j^J dxje-y^e^^y^ (4) 

where 

~ ' ^"\,.{T)-H{4>*Ar)AAr))) , (5) 



dr 13 



/3 = 1/T (T is the temperature), yj = ixj for the 25"- model or yj — for the 
4>S'-model, fi is the chemical potential and s is the state number per site allowed 
in each model. 

Now, it follows the standard procedures in which the replica method [5], 

Pn = --{\nZ{y})j^^ = Jnn^ (6) 

is used to get the configurational averaged thermodynamic potential. The repli- 
cated partition function {Z {y}"') j.^ is then given by: 

/oo POO 
n^9"7 / l[dpc.exp[Nn^{q^^,p^)] (7) 
-OO „, ^, — oo 



where TV = (/?J^iV/27r)"("+i)/2 anda = 1,2, - • • , n is the replica index. Within 
the static approximation, one has: 

(a, 7) a 

where the Fourier representation is used to express: 

^{2^> = n^/'"^^"<5~'"/ ^[Oa]exp[i?e//], (9) 

H^ff = J2 ^or + I 5] Pa^^ 5^ + 2 ^ q^^S^^S^ (10) 

a \ a (a,7) J 

with the convention that (0,7) indicates a distinct pair of replicas and the 
definitions: 

(p|^(a;) = ^0|^(a;) (/)|^(ijj)^ , = ±7r, ±37r, • • • . a;^ and are the Pauh matri- 



ces. 
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Therefore, in the IS-RSB Parisi scherae ^B], the thermodynamic potential 
is written from Eq. ([6]) as 



TO 

where 



- y Dzlnjy i:>i;[if(z,z;)]'"| -ln2 



(12) 



(z, i;) = cosh(/?M) + J DC cosh[v/A(z,t;,e)] , (13) 

with A(0,w,O = + (/3r)2 and 

Hz,v,£,) = jV2{y/q^z + y/qi - qav + ^/p - qiC), (14) 

where Dx = c?xe~^^/^/\/27r {x = z, v or ^). In Eqs. and (fH|) . go and gi 
are the IS-RSB order parameters and p = qaa = {Sa^a) is the diagonal repHca 
spin-spin correlation that is related to the diagonal replica matrix elements. 
The parameters go, ^i, P and to are given by the extreme condition of the 
thermodynamic potential Eq. p2p . The RS solution is recovered when go = 
qi{= q) and to = 0. 

Particularly, in the 4S model, the chemical potential /i controls the average 
occupation of fcrmions per site n = (n| + ni) which is obtained as: 

n ^ 1 + tanh(/3,) | (1 - p) - ipvf ( Dzl^^^^^'' ^ ^^t^^^' 



/ Dv[K{z,v)]"' 

(15) 

with 

cosh y/A{z,v,£,) sinh ^A(z,t;,,^) 
0(2;, C) = T7 TT -— • (lb) 

A{z,v,C) A^/^{z,v,Q 

The local stability of the RS solution for the FISG model is studied following 
close to de Almeida- Thouless analysis [24]. For this purpose, the stationary 
points (Eq. ([S])) are arbitrarily perturbed by independent quantities {S.a , rjaj) ■ 
The deviation AniS.a^Va-y) = ^niPa^qa-^) - ^n{Pa + , + Vai) cau bc 
obtained expanding A„ up to the second order in (^q,??^^): 

(17) 

a, 7 0,(71/) (q/3),(7I') 



where 

Go,7 = ^ - 16(5^^^5^5^) -I- 16(5^^^) (5^5^) = (A - B)5^,^ + B (18) 



Ga,7^ — -16(5^S'^S'^S'^) + 16(S'„5^)(S'^S'^) — <j b {i a ^'^-f and v ^^'^'^ 
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G 



2/32 J2 



P if a/3 = 



R. 



(20) 

The eigenvalues associated with the Hessian matrix (IT71) are known 26]. In 
the limit n ^ 0, there are three distinct eigenvalues: 



and 

where 



P-AQ + 3R = 



X±^[A- B + {P-4:Q + 3R)± Vu]/2 

Xat^P-2Q + R 

U=[{A-B)-{P-4Q + 3i?)]2 - 8(C - Df, 
1 



2/32 J2 



- / Dzyiiz)-4ip2{zM{z)+3ipt{z)], 



P-2Q + R = 
A-B^ 



1 



2/32 J2 
1 



- / Dz[^2{z)-vl{z)\ 



(21) 

(22) 

(23) 

(24) 

(25) 
(26) 
(27) 

(28) 
(29) 



In that case, the stable RS solution would occur when the Hessian eigenvalues 
A± and Xat are all positive. As one can see, different from the classical SK 
model, in the FISG model, the replica diagonal elements have an important role 
in the stability of the RS solution. Particularly, the condition for all eigenvalues 
of the Hessian matrix to be nonnegative in the paramagnetic solution (q = 0) 
is used to locate the tricritical temperature Ttc{T) at (for detail, see Ref. [14]): 



with 



- j Dz [fi{z) - (fliz)] , 
C-D = I Dz[ip2{z)^l{z)^^^{z)^^{z)\ 
d 



^C^(coshVA(z,C) 



/K{z), 



A{z, e) = {z, C) + /32r2 with h{z, = f3JV2i^z + Vp^O, and 
K{z) = ^—^ cosh/3^ + / coshJA{z,^). 



Ttc/J = -^^4/^2 
where ip2 and (^4 are defined in Eq. ()28p when q = and p = Ttc/JV2. 



(30) 
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3 Results 

The numerical solutions of the saddle point order parameters equations allow 
one to build phase diagrams (T/J) {T is the temperature) and F/J (see, for 
instance, Refs. [TH dS]). In the 45* model, it is also taken into account fj,/ J 
to build phase diagrams (see, for instance, Ref. [II]). The parameter J is 
related to the variance of the Gaussian random coupling given in Eq. 
Particularly, for the 45* case, the first order boundary has been located using the 
procedure introduced in Ref. [25] for the GS model. The Hessian eigenvalues 
(Eqs. (|21j|22p ) are inserted inside these phase diagrams. 

1 
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Figure 1: Phase diagram T/J versus F/J for the 2S model showing the behavior 
of the eigenvalues A+ and A_. The notation (+) ((— )) indicates a region where 
both eigenvalues A+ and A_ are real positive (negative). In the region indicated 
by i? > (ii < 0), these eigenvalues become complex with their real part 
positive (negative). The solid line indicates a second order transition. 

Fig. [1] shows a phase diagram T/J versus T/J for the 2S model. This is 
equivalent to study the problem in the spin space, as discussed previously. The 
SG phase boundary is followed by Almeida-Thouless line (defined as Xat = 0) 
which means that Xat is negative in the whole SG phase |15j. For very small 
T/J, eigenvalues A+ and A_ (denoted by A±) are real positive ( indicated by 
(+) in Fig. [1]). When F is enhanced, the contribution from the replica diago- 
nal matrix elements for the SG solution starts to increase. As a consequence, 
eigenvalues A± become complex conjugated pairs. In regions indicated by R 
(see Fig. [Ij, A± assume complex values with their real part positive {R > 0) 
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or negative {R < 0). At very low temperature, X± are real negative (indicate 
by (— ) in Fig. [T]). The onset of regions with R > and R < always occurs 
below Xat — 0. However, when T increases, boundaries of such regions, as well 
as that one with (— ), become increasingly close to the AT line. In fact, these 
boundaries go toward the QCP given by Fc = 2 J, which suggest that the posi- 
tivity of Xat in the FISG model could not be enough to guarantee the replica 
symmetry stability near the QCP. 




Figure 2: Phase diagrams T/J versus T / J showing the behavior of eigenvalues 
A+ and A_ for the AS model and several values of /x/J. It is used the same 
convention as Fig. [1] for the instability regions of the RS solution. The solid 
lines indicate second order phase transition. Dashed lines indicate a first order 
phase transition, in which results within RS and IS-RSB solution are compared. 
Tfc indicates the tricritical point. 

In the 45 model, the combined effects of /i and F on Xat and the X± can 
be seen in Fig. [51 This figure shows phase diagrams T/J versus T/J for sev- 
eral values of /i/J. For all values chosen for /i/J, the freezing temperature Tf 
decreases toward a QCP when F — > F^, as in the 25-model. The half-filling sit- 
uation = 0) is exhibited in Fig. [2la). In this case, the SG solution is always 
unstable against the replicon stability analysis. However, in contrast with the 
2S-model, it now appear a region with R > even for F = 0. Fig. 2(b) presents 
the same characteristics found in Fig. 2(a), except there is no region with (+) 
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Figure 3: Phase diagrams T/J versus [ij J for the 4S model with for several 
values of r/J. The same convention as Fig. [T]is used for the instability regions 
of the RS solution. The solid lines indicate second order transition while dashed 
lines indicate a first order transition, in which the RS and IS-RSB solutions 
are compared. 



at low T/J and T/J. Nevertheless, for higher values of /i/J, as shown in Figs. 
2(c)-(d), the consequences for SG instability from the interplay between charge 
and spin correlations become more important. For instance, when fi/ J ~ 1 (see 
Fig. d^c)), there are no SG solutions for < F < 0.5/ J. When F > 0.5J, 
with the arising of these solutions, there is also a change in the nature of the 
SG boundary phase with the onset of a tricritical point at temperature Ttc. 
Below Ttc, a first order boundary is located in RS and IS-RSB schemes, but, 
without any important difference between them [14]. Regarding the behavior 
of Xat at high values of fi/ J, it is found that the replicon remains negative in 
the entire SG phase, as previously observed in Figs. 2(a)-(b). However, there 
is an enlargement of regions with R > 0, R < and (— ). Particularly, close 
to the first order boundary phase, there are only regions with i? < or (— ). 
On the contrary, close to the second order boundary phase, the region with 
(+) becomes smaller. Finally, in Fig. [HJd), the SG phase boundary is strongly 
decreased as compared with Fig. ^c). The important point in this last case is 
that the regions with i? < or (— ) occupy almost the whole SG phase. 
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Figure 4: Average occupation n versus fi/J for several values oiT/J and T/J 
of the SG solution. The dashed-dotted lines indicate the IS-RSB solution. The 
insets show in detail the difference between the value of n in RS and IS-RSB 
solution. 



The combined effects of F and fi can also be seen in a different perspective. 
In Fig. [3l the phase diagrams are built now in T/J-fi/J surfaces for several 
values of F/J. Fig. [3^a) shows a phase diagram for F with regions R > 0, 
R < and (— ). This result can be compared with Ref. [25j using the mapping 
introduced in Ref. [13j . The phase boundary transition at T/ is continuous at 
low fi/J. When /i/J increases, a tricritical point appears at Ttc and, at lower 
temperature, the first order part of Tf presents a reentrance (see discussion in 
Ref. [Mj). Close to first order boundary, besides the usual replicon instability, 
the eigenvalues A± become complex with i? < as shown in Fig. [2l The 
enhancement of F, shown in Figs. [3I^b)-(d), decreases the second order boundary 
phase and the tricritical point as well as it suppresses the reentrance in the first 
order one [M]. The behavior of Xat remains basically the same found in Fig. [2l 
The SG phase is always replica symmetry breaking. However, again, the region 
with (+) decreases, while those ones with R > 0, R < and (— ) increase. It 
should be emphasized that, close to the first order transition, it is observed, 
again, a region where, at same time, the replicon and the X± violate the replica 
stability condition. 

The complicated connection between charge and spin correlations due to the 
presence of F can also be seen in the average occupation of fermions per site n 
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(see Eqs. (HHl-dlll)). Figure H shows the behavior of n within RS and IS-RSB 
schemes as a function of /i/J for several isotherms and F/J. For F = 0, when 
temperature is decreased, n remains at the half-fiUing even for fi/ J ^ 0. This is 
consistent with an earher result [ini HT] , which suggests that the ocuppation of 
nonmagnetic states is exponentially small when T is decreased. However, it is 
also true that this effect is weakened when IS-RSB solution is used, particularly, 
at lower T [13]. Nevertheless, the presence of F enforces the tendency to preserve 
the half- filling occupation, as can be seen, for higher T in Fig. |DJa). Indeed, 
this effect is enhanced when temperature is lowered as shown in Figs. [IJb)-(d). 
Most important, when F is turned on, the differences between RS and IS-RSB 
solutions become irrelevant. In other words, F enforces the effect of temperature 
preserving n at the half-filling for both RS and IS-RSB levels of description. 

4 Conclusions 

In this work, it has been studied the stability of spin-glass phase in the FISG 
model when a magnetic transverse field F is applied and, particularly, X± (the 
longitudinal eigenvalues of the Hessian matrix) . The reason is that the presence 
of F produces important effects in the diagonal replica matrix elements which 
are closely related to the behavior of \±. In fact, the FISG model can be 
presented in two formulations, the so-called AS and 25* models [TB]. In the 
45* model, the original four states of the Fock space are maintained. In the 2S 
one, the nonmagnetic states are not allowed. Particularly, this means that in 
the 45' model there is a connection between charge and spin degrees of freedom 

The results presented in Figs. [TJ2] display two distinct situations for Xat 
and A±. For the 2S and 45 models, the results show that the replicon Xat is 
always real negative in the entire SG phase for any value of F/J. On the other 
hand, for the 2S model (see Fig. [T]), the enhancement of F/J produces regions 
where A± becomes complex with real part positive (i? > 0), negative {R < 0) 
or even real negative ((—)). Furthermore, all boundary lines of such regions 
converge to the QCP when F ^ Fc. In the 45 model (see Fig. [5]), X± also 
present regions with R > 0, R < and (— ). In fact, there are other similarities 
between the two models. When F — > Fc, the boundary lines of these regions 
have the same behavior already found in the 25 model. Most important, this 
result is independent of a particular value of ^/J. Nevertheless, for F < F^ 
when < /i < J, the regions with R < and (— ) are relatively small and they 
appear only at lower temperature. By contrast, the situation is different when 
/z > J. Particularly, regions with R < and (— ) are considerably enlarged while 
that one with (-|~) is decreased. It should be remarked that close to the first 
order boundary phase, besides the region with (— ), there is only region with 
R<0. 

The role of F described in the previous paragraph for the 45 model is con- 
firmed in Fig. [3l There, the increase of fi produces regions with R > 0, R < 
and (— ) as expected, for instance, from the mapping with the GS model when 



12 



r = [13]. However, the increase of T tends to depress Tf, the tricritical point 
Ttc and to destroy the reentrance in the first boundary phase [Mj . Such increase 
also depresses boundary hues of regions with R > 0, R < and (— ). The im- 
portant point is that these boundary hnes follow the change in location of the 
Ttc- Therefore, no matter the value of /x, the increase of F tends to reproduce 
the scenario in which all boundary lines of regions R > 0, R < and (— ) are 
decreasing towards zero temperature. 

The behavior of n exhibited in Fig. 2] can clarify the role of F in the problem 
as source of similarities between the 25' and 43 model when F Tc- The 
increase of F leads n to remain at the half-filling even when /x is increased from 
zero. Therefore, the increase of F acts to avoid the double occupation. In the 
opposite limit, when F << Fc, the increase of chemical pressure tends to remove 
n from the half-filling which leads the 45 model to become closer to the classical 
GS model. 

Therefore, it is possible to identify two regimes. In the first one, which is 
common to both 25 and 45 models, the region with (+) is increasingly small 
when F — > Fc. In the second regime, when F is not close to Fc, there is a mixing 
of effects due to the presence of fi and F. In that case, the increase of /x produces 
the increase of regions where the eigenvalues X± are complex or negative. It 
should be remarked that the regions with R < and (-) are exactly those ones 
which occupy most part of the SG phase diagram, particularly, close to the first 
order boundary phase. Therefore, even as suggested in Refs. [121 [3D], that the 
condition R > would be enough to garantee the SG stability, the dominance 
of regions with R < and ) for certains values of F and /i indicate clearly 
an additional instability to the SG solution which must be taken into account 
simultaneously with that one related with the replicon. 

Surely, the reliability of some results are quite limited by the use of SA, 
particularly, those ones close to the QCP. Nevertheless, in the first regime, the 
results show clearly that the tendency to decrease the region with when 
F is increased, can be observed even at higher temperature. This tendency 
would indicate that close to the QCP the analysis of X± could also play a more 
important role. However, there is no doubt that it is needed to go beyond the 
SA to obtain a more conclusive result close to the QCP. On the other hand, for 
the second regime, results are more conclusive. The effects of F combined with 
the increase of /i give an additional instability at higher temperatures. This new 
instability is important, particularly, close to the first order boundary where the 
use of SA is not so limited as in the case of a continous transition [T9^. 

To conclude, in this work, it has been studied in detail the behavior of X± 
in the FISG model as a function of F. It has been presented several evidences 
indicating that these particular eigenvalues of the Hessian matrix can be also 
source of problem to the local stability of the RS solution when F increases. It 
is clear that this results are very much restricted to the FISG model. However, 
one can speculate whether or not others diferent formulations of infinite-range 
Ising spin-glass model in a transverse magnetic field could experiment similar 
additional instabilities likewise those ones found in the 25 model. 

One last comment must be done. Although there are strong similarities 
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between the 25 and 45 models in the first regime, it is not obvious that this 
similarity can be explained by the behavior of the n which remains at the half- 
filling when r increases. One alternative, would consist of adding to the model 
given in Eq. ([1]) a local on site repulsive Coulomb interaction with strength U. 
Therefore, in the limit U oo, it would be possible to avoid locally the double 
occupation [13]. This approach is currently under investigation. 
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